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CO adsorption on an Au monatomic chain is studied within density functional theory in nanocon- 
tact geometries as a function of the contact stretching. We compare the bridge and atop adsorption 
sites of CO, finding that the bridge site is energetically favored at all strains studied here. Atop 
adsorption gives rise to an almost complete suppression of the ballistic conductance of the nanocon- 
tact, while adsorption at the bridge site results in a conductance value close to 0.6 Go, in agreement 
with previous experimental data. We show that only the bridge site can qualitatively account for 
the evolution of the conductance as a function of the contact stretching observed in the experimental 
conductance traces. The numerical discrepancy between the theoretical and experimental conduc- 
tance slopes is rationalized through a simple model for the elastic response of the metallic leads. We 
also verify that our conductance values are not affected by the specific choice of the nanocontact 
geometry by comparing two different atomistic models for the tips. 



PACS numbers: 81.07.Lk, 73.63.Rt, 73. 23. Ad, 73.20.Hb 



I. INTRODUCTION 

As demonstrated by a recent experiment^ the atomic 
conductance of a Au nanocontact can be altered by ad- 
mitting gaseous CO in proximity of the nanocontact, sug- 
gesting that CO chemically interacts with the Au atoms 
in the thinnest part of the contact (a monatomic chain 
or a single atom). Indeed, besides a slight renormaliza- 
tion of the typical conductance peak at 1 Go, the conduc- 
tance histogram of Au after exposition to CO displays 
an additional peak between 0.5 Go and 0.6 Go, which can 
be attributed to CO adsorption. This is further sup- 
ported by the analysis of the conductance traces, which 
clearly pinpoints the existence of an atomic structure 
with a CO molecule adsorbed on, or incorporated into 
the monatomic chain. This structure is invariably char- 
acterized by a conductance value close to that of the new 
histogram peak and also by a small increase of the con- 
ductance upon further pulling of the nanocontact, just 
prior to contact breaking^ 

In a previous density functional study, we addressed 
the adsorption of CO on Au monatomic chains using an 
infinite chain geometry without tips. 2 The bridge adsorp- 
tion site was found to be energetically favored with re- 
spect to the atop site, both at the equilibrium spacing 
of the Au chain and at larger values of the Au-Au spac- 
ing. We characterized the adsorption process by identify- 
ing the bonding/antibonding pairs of 5a and 2ir* states 
which arise from the hybridization between CO molecu- 
lar levels and Au metal states. Therefore, the electronic 
structure of this chain/adsorbate system results from a 
donation/backdonation mechanism analogous to the Bly- 
holder model for CO on transition metal surfaces^ - — as 
previously demonstrated also for CO on Pt monatomic 
chains^ 7 - 



A strong connection between some electronic structure 
features of the CO/Au-chain system and the effects pro- 
duced by CO adsorption on the ballistic transport across 
the chain was established in that work.— In particular, the 
coupling of the 5a antibonding state with the s states of 
the Au chain was found to generate a dip in the s trans- 
mission. Hence, the position of this state turns out to 
be crucial in determining the reduction of the (tipless) 
conductance due to the impurity, as well as the strain 
dependence of the conductance in presence of CO. In 
the atop geometry, the 5a antibonding state is located 
very close to the Fermi level (Ep) and approaches Ep 
as the Au chain gets stretched, hence it strongly reduces 
the conductance of a moderately strained chain and cuts 
almost completely the conductance of a highly strained 
chain (close to the rupture Au-Au spacing) <2&2 In the 
bridge geometry instead, the conductance reduction is 
much lower than in the atop geometry because the 5a 
antibonding state is more distant from Ep at medium 
strains and the associated transmission dip actually dis- 
appears at high strains 3^ At variance with the atop 
geometry, which shows a decrease of the conductance 
with strain and is therefore not compatible with the ex- 
perimental conductance traces, the conductance of the 
bridge geometry shows the correct dependence on strain, 
namely, a slight increase with increasing strain^ 

The calculated conductances of the tipless bridge ge- 
ometry are however too large compared to the experimen- 
tal conductance of the histogram peak appearing upon 
CO adsorption.i Moreover, the nonuniform distribution 
of the strain, the finite length of the chain, and the ef- 
fect of the tips are not described in the infinite chain 
model, therefore those results need to be corroborated 
by a more realistic modeling of the nanocontact. Density 
functional calculations for various kinds of impurities ad- 




FIG. 1. (Color online) Lateral view of the periodic cells for 
the Au chain between Au(OOl) surfaces and CO adsorbed at 
the (a) bridge and (b) atop sites. The smaller cell in (a) is 
used to compute the complex band structure of the leads. 
The spacing between the bulk Au layers (db), the outermost 
interlayer distance (d s ), and the inter-surface distance (d aa ) 
are indicated. 



sorbed on model Au nanocontacts (for instance, molecu- 
lar hydroge n 11 ' 12 and oxygen^), for CO on other metals, 
such as Ptr^ and also for one specific geometry of CO on 
a Au nanowire contact^ are already available in the lit- 
erature. However, different adsorption geometries of CO 
on Au nanocontacts have not been compared and the ef- 
fect of contact stretching on the ballistic conductance of 
this system still needs to be analyzed. 

In this work, we study the adsorption energetics of CO 
on a short Au monatomic chain and we examine the bal- 
listic conductance of the system as a function of the con- 
tact stretching using model nanocontact geometries. By 
comparing the bridge and atop adsorption geometries of 
CO, we find that the bridge site is energetically favored at 
all strains, confirming the result obtained within the infi- 
nite chain model.- In the atop geometry, the Au conduc- 
tance cut due to the ha antibonding states^ is observed 
also in the present nanocontact model when the contact 
is sufficiently stretched. The bridge geometry instead, 
reproduces to a good accuracy the experimental conduc- 
tance value (0.5 -j- 0.6 Go)r representing a significant im- 
provement with respect to the infinite chain model. The 
slope of the conductance as a function of the contact 
stretching is instead reproduced only qualitatively by our 
nanocontact model. We address this discrepancy show- 
ing that a conductance slope closer to the experimental 
one can be obtained by taking into account, even only 
approximately, the effect of the elastic response of the 
bulk leads to the external pulling force, not described by 
our nanocontact model. 

This paper is organized as follows: in Sec. |nl we 
describe the numerical methods and approximations 
adopted in the calculations; in Sec. IIIII and in Sec. IIVI 
we will report on, respectively, the adsorption energetics 
and the ballistic conductance of the system as a function 
of the contact stretching; Sec. [V] is devoted to a discus- 



sion of the conductance results in the light of the available 
experimental data; finally, our conclusions will follow in 
Sec. ED 



METHODS AND COMPUTATIONAL 
DETAILS 



II. 



All density functional calculations presented in this 
work are carried out using the plane-wave pseudopoten- 
tial code PWscf contained in the Quantum ESPRESSO 
packaged The exchange and correlations functionals, 
both for the local density (LDA) and the generalized 
gradient approximation (GGA), the basis set cutoffs, the 
pseudopotentials, and the smearing parameters are the 
same as those used for the infinite chain model^. 

The monatomic chain in the nanocontact geometries 
is modeled as a row of metal atoms suspended between 
fee bulk leads terminated by two facing (001) surfaces. 17 
In Fig. [I] we show the simulation cells consisting of a 
periodically repeated slab geometry with 7 fcc-Au layers 
perpendicular to the [001] direction, which is chosen as 
z axis and coincides with the electron transport direc- 
tion (see later). The spacing between the inner layers 
corresponds to the theoretical equilibrium value in the 
bulk d b , which is 2.029 A (2.082 A) according to our LDA 
(GGA) calculations. The clean unreconstructed Au(001) 
surface shows a significant inward relaxation of about 
1.8% (1.5%) for outermost layer within the LDA (GGA), 
in agreement with recent DFT calculations, 18 while spac- 
ings between inner layers deviate from the bulk value by 
less than 0.4%. Therefore, we keep into account the re- 
laxation of the outermost layer on both sides of the junc- 
tion by setting the first interlayer distance d s = 1.992 A 
(2.051 A). The in-plane distances between gold atoms are 
those of the bulk. The short chain is freely suspended in 
the vacuum region between the slabs and the apex atoms 
of the chain are attached to the surfaces at a 4-fold coor- 
dinated hollow site. The symmetry group of the straight 
chain between the two (001) surfaces is D^, the same 
as in the clean surface case studied with a slab geometry 
having inversion symmetry. 

A single CO molecule is adsorbed at the center of 
the Au chain, at the bridge site of a 4-atom-long chain 
(Fig. QJi) and at the atop site of a 5-atom-long chain 
(Fig. QJ>). The symmetry group of the system in pres- 
ence of CO is C2v, as in the case of CO adsorbed at 
the bridge or atop sites of the infinite chain. The xy in- 
plane periodicity of the simulation cells corresponds to a 
(2\/2x 2-\/2)R45 surface structure, which gives a chain- 
chain spacing of about 8.12 A in the x and y directions 
between two adjacent replicas. The full Brillouin zone 
(BZ) of these structures is sampled with a uniform mesh 
of 6 x 6 x 3 k-points, which can be reduced by symme- 
try to 12 points (18 when CO is adsorbed). We relax 
the atomic positions of the Au chain and of CO until the 
forces on those atoms drop below 0.026 eV/A. We also 
check that the optimized atomic positions of the chain 



and of CO obtained in this way result in atomic forces 
below 0.08 eV/A when used in a (3a/2 x 3^2)1145° cell 
(chain-chain spacing of 12.18 A along x and y). For this 
larger cell we reduce the size of the k-point mesh down 
to 4 x 4 x 3. 

The ballistic conductance is evaluated with the 
Landauer-Buttiker formula, G = e 2 /h T(Ep), where 
T(Ep) is the total transmission at the Fermi energy. We 
calculate the electr on transmission using the scattering- 
based approach of IChoi and IhmU ^ extended to ultra- 
soft pseudopotentials^ and implemented in the PWcond 
code. 10 In transmission calculations, we include in the 
scattering region a portion of the leads together with 
the impurity region (see Fig. []}. The semi- infinite bulk 
leads are modeled using an additional, smaller cell which 
coincides with two bulk Au(001) layers (smaller cell in 
Fig.CQi). 21 In the (2a/2 x 2\/2)R45 nanocontact geome- 
try (which we will also call "abrupt" junction), the k^- 
dependent transmission is sampled with a uniform 7x7 
shifted mesh of k^-points in the 2D-BZ perpendicular to 
the transport direction, corresponding to 10 and 16 k^- 
points in the irreducible 2D BZs of the clean nanocon- 
tact and of the nanocontact with CO, respectively^ 
Smoother junctions have been simulated using periodic 
cells with a (3 x 3) in-plane periodicity and pyramidal 
tips connecting the chain apexes to the Au(001) surfaces 
(see Fig. [4]). This geometry is composed of a seven-layer 
slab with an interlayer spacing equal to dh, plus four ad- 
ditional Au atoms in the positions of an additional layer 
at a distance d s from the surface planes, with the chain 
attached to the 4-fold hollow sites formed by these ad- 
ditional atoms. A 5 x 5 uniform mesh of k^-points has 
been used to sample the transmission in the 2D-BZ. 



TABLE I. Optimized distances (in A) and chemisorption en- 
ergies (in eV) of CO at the bridge site of a 4-atom-long chain 
(cf. Fig. [T^) obtained within GGA for selected values of d ss . 
The optimized distances of the clean chain are also reported. 
Owing to the symmetry, a?3_4 = di_2 and d4 S = d s i. 



III. GEOMETRY AND ENERGETICS 



For different values of the inter-surface distance d ss (cf. 
Fig. [p, we consider straight monatomic chains without 
CO and chains with CO adsorbed at the bridge or atop 
sites. In this way, by increasing d ss we mimic the in- 
crease of strain produced on the chain by the pulling of 
the contact ends. We optimize the atomic positions of 
the C and O atoms, and of the Au atoms belonging to 
the chain, while the positions of all Au atoms in the (001) 
planes are kept fixed (as are the distances between the 
planes) . Since we are mainly interested in the local inter- 
action between CO and the Au chain, which is primarily 
affected by the chain strain, we will not consider here 
the atomic relaxations of the Au planes at the two sides 
of the junction. The effects of these relaxations on the 
strain dependence of the ballistic conductance are stud- 
ied in Sec. fVl through an approximate model of the elastic 
response of the atomic planes to the contact stretching. 





4-Au chain 




CO at the bridge site 




dss 


dsl dl-2 ^2-3 


dsi 


dl-2 Gfo-3 dAu-C dc-o 


-C'chcm 


11.56 


1.87 2.61 2.61 


1.90 


2.59 2.86 2.03 1.17 


-1.47 


12.16 


1.99 2.72 2.74 


1.93 


2.61 3.09 2.02 1.17 


-1.61 


12.76 


2.09 2.82 2.94 


1.99 


2.67 3.44 2.06 1.18 


-1.94 



A. Clean nanocontact geometries 

We first consider a straight chain without CO and op- 
timize the z coordinate of the Na u Au atoms of the 
chain, thus removing the constraint of uniform inter- 
atomic spacing that was adopted in the infinite chain 
model.— We did not explore here zigzag or bent config- 
urations, which are expected to become favored only at 
low values of d ss in the clean cham^r— The optimized 
distances between the Au atoms in the chain (di-2, ^2-3, 
. . . ) and the distance between the surface plane and the 
apex atom of the chain (d s i) are presented in the left 
parts of Table [J and Table [TT| for a 4-atom-long and a 5- 
atom-long chain, respectively. We report here only GGA 
data, since the LDA results give the same qualitative 
pictured In the first row of each table, d ss is chosen 
to give Au-Au distances in the short chain close to the 
equilibrium value in the infinite chain (2.61 A), while in 
the second and third rows the selected d ss values result 
in moderately or highly stretched Au-Au bonds, respec- 
tively. 

At the lowest strain considered here, the atoms in the 
chain are almost equally spaced, but the Au-Au bond 
length at the extremities of the chain, g?i_2, adjusts to 
a value slightly smaller than in the middle (only c?2-3 in 
the 4-atom and 5-atom chains, but we verified that this 
holds also for the inner bonds of longer chains) £L As the 
nanocontact is stretched, the distance di_2 becomes pro- 
gressively shorter than ^2-3- By studying longer chains 
(Aau = 6 and Aa u = 7, not reported here), we observe 
Au-Au bond lengths which increase while going from the 
ends toward the center of the unstrained chains. Instead, 
when the average Au-Au bond length is above 3.0 A the 
chains show a tendency to dimerization, with alternat- 
ing longer and shorter bonds, as already reported in the 
literature.— 



B. Geometries with adsorbed CO 

We now consider the nanocontact geometries with an 
adsorbed CO molecule. The atomic structure is partially 
optimized as described above and the positions of the Au 
atoms in the chain are fully relaxed. We first consider the 



TABLE II. Optimized distances (in A) and chemisorption en- 
ergies (in eV) of CO at the atop site of a 5-atom long chain 
(Fig. [Tp). Here, d 4 _ 5 = di_ 2 , d 3 -4 = cfe-3, and d 5s = d sl . 





5-Au chain 




CO at the atop site 




Cf'ss 


d s l di-2 ^2-3 


dsl 


rfi-2 d'2-3 <2au-c dc-o 


-L'chcm 


14.16 


1.86 2.61 2.61 


1.98 


2.72 2.76 1.98 1.14 


-0.78 


14.96 


2.00 2.73 2.75 


1.99 


2.73 2.83 1.99 1.14 


-0.55 


15.76 


2.09 2.87 2.92 


2.04 


2.77 3.08 1.98 1.14 


-0.63 



wards the molecule and the two lateral atoms, Au(2) and 
Au( 4 ), slightly displaced downwards to create a zigzag 
geometry. By comparing di_2 with d<i— 3, we see that 
the binding between C and Arj.(3) weakens the metallic 
bond between Au(3) and the two neighbouring atoms, 
Au( 2 ) and Au(4). Indeed, when CO is adsorbed atop the 
ratio d2—s/di-2 is always larger than in the pristine 5- 
atom-long chain and grows more rapidly with strain (cf. 
Table HJ). 



bridge geometry, where the CO is placed upright at the 
bridge site between the two central Au atoms (Au( 2 ) and 
Au( 3 )), as shown in Fig. [T^. In Table fl] (right side), we 
report the optimized distances for the three d ss values 
considered before. The carbon-oxygen bond length in 
the adsorbed molecule (dc-o) is equal or slightly larger 
than in the infinite chain model 2 (by less than 1%), while 
C-Au bond lengths (cJau-c) are 1% to 3% larger. We also 
notice that the bond length between the two Au atoms in 
contact with CO ((^2-3) is always longer than the other 
Au-Au bonds in the chain (di_2 = (^3-4). This Au-Au 
bond softening in correspondence of the adsorption site 
is observed also for longer chains^ At the lowest strain 
studied here (d ss — 11.56 A), d 2 -3 adjusts to a value 
similar to or slightly larger than the Au-Au distance cor- 
responding to the bridge energy minimum for the infinite 
chain geometry 2 - at g?au-Au = 2.87 A. The other Au-Au 
bond lenghts, instead, are closer to the equilibrium spac- 
ing of the infinite chain. At low and moderate strains, the 
chain in its relaxed geometry bends towards CO, while 
at larger strains the chain atoms get progressively more 
aligned forming an almost linear strand with one over- 
stretched Au-Au bond in correspondence of the adsorp- 
tion site. Because of the flexibility of the Au-C-Au bond 
angle, very large Au-Au distances between the two atoms 
in contact with CO are possible as d ss is increased, but 
the tilting of the molecule axis from the perpendicular 
position could become favorable above a critical value of 
d ss ^ However, the upright bridge position in the infi- 
nite chain geometry 2 - is the lowest energy configuration 
for Au-Au distances up to 4.2 A, a distance longer than 
the largest ^2-3 reported in Table |U Hence, one can ex- 
pect that the perpendicular position of CO is preferred 
for all values of d ss considered here. 

In Table HIl we report the optimized distances for CO 
adsorbed atop the central atom of a five-atom-long Au 
chain (Au( 3 \ in Fig. [Tp). The carbon-oxygen bond dis- 
tance dc-o changes very little with strain and is very 
close to the value found in the infinite chain model. 2 Also 
the C-Au( 3 ) bond length is not much influenced by strain 
and is only about 1% larger than the corresponding dis- 
tance for CO on the infinite chain. As can be seen from 
the two nanocontact geometries reported in Fig. [TJ the 
adsorption of CO in the atop position gives rise to larger 
distortions of the Au chain compared to the bridge posi- 
tion, with the atom which binds to C, Au(3), moving to- 



C. Chemisorption energies 

We now turn to examine the chemisorption ener- 
gies of CO as a function of the surface-surface distance 
E c hem(d ss ), which are reported in Table [J and Table [TT1 
for the bridge and atop geometries, respectively. The 
large difference between the chemisorption energies of the 
bridge and atop adsorption sites indicates a strong pref- 
erence for the former, as already observed in the infinite 
chain models In the bridge geometry, the chemisorption 
energy decreases with strain, following the same trend 
seen for the infinite chain model; 2 - while in the atop ge- 
ometry this happens only at large enough values of d ss . 

This can be seen more clearly in Fig. [21 where we 
compare the chemisorption energies computed with the 
nanocontact geometry as a function of d ss , -E c hcm(<4s), 
with those obtained from the infinite chain geometry as 
a function of the uniform spacing of the chain d c hain, 

Schom(rfchain)^ To do this, we express -Echcm(dchain) 

as a function of an equivalent inter-surface distance 
d S s(d c ha,in) , which can be directly compared to d ss . Hence, 
we define: d ss = 2 • (d s i} + (Na u — 1) • ^chain, where the Au- 
Au spacing <i c hain is multiplied by the number of Au-Au 
bonds in the A^u-atom-long suspended chain and (d s i) 
accounts for the distance between the surface plane and 
the apex atom of the chain. We obtain (d s i) by aver- 
aging the d s i values reported in Table fl] for the bridge 
geometry and in Table [ill for the atop geometry, which 
give (d s i) = 1-94 A and (d s i) = 2.00 A, respectively. 

For the bridge geometry, there is a qualitative agree- 
ment between -E c hcm(<is) from the infinite straight chain 
(open squares in Fig. [2]) and E c h cnl (d ss ) from the 
nanocontact geometry with Na u = 4 (filled squares). 
This difference in the numerical values could be mainly 
imputed to the nonuniform Au-Au spacing, to the bend- 
ing, or to the short length of the chain in the nanocontact 
geometry. For the atop geometry (open and filled cir- 
cles), there is a larger deviation at low d ss because of the 
appearance of zigzag configurations in the nanocontact. 
When the height of the Au atom below CO is optimized 
also in the infinite chain (open diamonds) , the agreement 
with the chemisorption energy obtained in the nanocon- 
tact improves considerably, because that Au atom shows 
a pronounced displacement toward the CO at low strains. 
Consequently, energy contributions to -E c hem due to dis- 
tortions of the chain further away from the adsorption 
site should be of smaller importance. 
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FIG. 2. (Color online) Chemisorption energies E c hem{d as ) for 
CO at the bridge site of a four-Au-atom chain (filled squares) 
and at the atop site of a five-Au-atom chain (filled circles). 
The values of -E c hem(d c hain) for the infinite straight chairs are 
reported here as a function of ri ss (ri c hain ) = 2-(d s i) + (A r Au — 1)- 
dchain, where JVa u = 4 for the bridge geometry (open squares) 
and JVau = 5 for the atop geometry (open circles), (dsi) is 
chosen as described in the text. For the atop site of the infinite 
chain we also report E c h cm (d ss ) after optimizing the position 
along x of the Au atom below CO (open diamonds). 



Another point that influences the numerical values of 
E c hem is the choice of the clean nanocontact geometry 
which gives the reference energy entering into the calcula- 
tion of -E c h cm . Zigzag and bent configurations have been 
predicted for short monatomic chains between tips, 23 
while these have not been considered in the present work. 
However, we expect that only the -E c hem values for the 
smallest d ss value considered here would be affected if 
using a nonlinear chain geometry as reference. Indeed, 
zigzag or bent geometries are preferred to the linear chain 
only for sufficiently small values of d ss , while the linear 
chain is favored with respect to both zigzag and bent 
configurations when the average spacing between atoms 
in the chain is equal or larger than about 2.65 A and 
the three configurations are nearly degenerate for slightly 
smaller spacings. 23 

Finally, this comparison between the infinite chain and 
the nanocontact geometry shows that a very simplified 
model such as the straight infinite chain with uniform 
spacing gives chemisorption energies of CO in qualita- 
tive agreement with those obtained through a more re- 
alistic and computationally expensive model, leading to 
the same site preference prediction. A more quantitative 
agreement in -E c hem can be obtained by including a few 
additional degrees of freedom in the geometry optimiza- 
tions of the infinite chain, such as the Au atom displace- 
ment when CO is atop, or a larger spacing between the 
two Au atoms in contact with CO in the bridge position, 
which could be inferred from the residual forces after a 
partial structural optimization. 



IV. BALLISTIC CONDUCTANCE 

In this section, we address the effects of CO adsorp- 
tion on the electron transport of the Au nanocontact for 
several values of the inter-surface distance, d ss . In this 
way, we aim to describe the strain dependence of the 
CO-induced changes in the ballistic transmission and to 
simulate the evolution of the conductance as the contact 
gets gradually stretched. 



A. Clean nanocontact 

We start discussing the ballistic transport of the short 
chain for the abrupt junction without CO. The ballis- 
tic transmission as a function of the electron scattering 
energy is presented here for the LDA case only. How- 
ever, the GGA ballistic conductances reported later in 
Sec. |V] (Fig. [5]) show that the same conclusions could be 
reached within GGA. In Figs.G3a-b), we show the trans- 
mission of a 4Au-atom straight chain as a function of the 
scattering energy for two selected values of the intersur- 
face distance d ss , one resulting in a moderately strained 
chain [d ss = 11.76 A, Fig.^a)], and the other in a highly 
strained chain [d ss = 12.36 A, Fig.^b)]^ 

In the low strain configuration [Fig. [3£a)], the conduc- 
tance is slightly above 0.9 Go, in fair agreement with the 
theoretical conductance of a 4-atom-chain between Au 
tips oriented along the [110] direction.— Although the 
number of scattering channels increases with the cross- 
section of the superccll, for a long enough chain and a 
small charge transfer between the leads and the chain, 
the theoretical maximum of the transmission is given by 
the number of channels in the infinite tipless chain. Thus, 
the conductance value close to 1 Go can be attributed to a 
single well transmitted spin-degenerate channel of s char- 
acter. Around the Fermi energy (Ep) the transmission 
curve is rather flat, while below Ep it is more structured 
because of a slightly higher reflection of the s channel 
and the additional contribution from the poorly trans- 
mitted d channels of the chain. This is in line with the 
common understanding that ri-type scattering states are 
more reflected by the presence of an abrupt change in 
the atomic geometry, like that at the surface/wire inter- 
face, owing to the highly directional character of the d 
wavcfunctions. 

When a larger strain is considered [Fig.[3[b)], the con- 
ductance changes very little and decreases of just about 
1 %. Previous calculation s 23 ' 27 have also found that the 
conductance decreases monotonically when the junction 
gets stretched and that the conductance has not yet 
started to drop significantly for average Au-Au spac- 
ings around 2.80 A, being still close to 1 Gor 7 . or slightly 
below i^ 3 . Above Ep the shape of the transmission func- 
tion is very similar to that obtained for the low-strain 
configuration, while below Ep the reflection slightly in- 
creases with strain. 

In Figs. 131 c-d), we show the transmission function for 
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FIG. 3. (Color online) Transmission as a function of the electron scattering energy for the clean Au nanocontact (diamonds) 
and for the Au nanocontact with CO (squares). On the left (right) panels, we report the transmission of the clean Au chain 
with A'au = 4 (JVau = 5) and of the chain with CO adsorbed at the bridge (atop) site. For each system, two values of the 
distance between the Au(OOl) surfaces are considered, a smaller one, corresponding to a moderately strained chain (top panels), 
and a larger one imposing a high strain on the chain (bottom panels). In each plot, the Fermi energy is indicated by a solid 
vertical line and the central part of the scattering region with the impurity is shown in the insets. The inter-surface distance 
d ss and the optimized distances (in A), as well as the chemisorption energies of CO (in eV) are also indicated. 



the 5-atom-long chain for d ss = 14.56 A and for d ss = 
15.36 A, corresponding to moderate and high strains, re- 
spectively. The conductance of the 5Au-atom chain is 
about 0.96 Go for both d ss values considered here and 
therefore is slightly larger than that of the 4Au-atom 
chain, in agreement with the odd-even effect seen in 
experiments 29 and in theoretical calculations ! 25 ' 30 For in- 
stance, de la Vega et al— have calculated the theoretical 
conductance of short chains between flat Au(lll) sur- 
faces showing that the conductance for an even number of 
atoms in the chain is lower than that for an odd number 
of atoms, as also seen in experiments. Their conductance 
values, about 0.96 Go and 0.99 Go for 4-atom and 5-atom 
long chains, respectively, are somewhat larger compared 
to the results of this work, probably because their chains 
are attached to the more compact Au(lll) surface. Also 
for the 5Au-atom chain, we find that the conductance 
dependence on strain is quite modest, since it stays al- 
most constant when the average spacing in the chain 
grows from about 2.74 A to 2.90 A as d ss is increased. 
For the 5-atom-long chain between Au(lll) surfaces, de 



la Vega and coworkers have found that the conductance 
goes from 0.99 Go to about 1 Go when the spacing be- 
tween atoms in the chain grows from 2.70 A to 3.00 A, 
while the conductance of the 4-atom-long chain does not 
change appreciably— as in the case of the Au(001) sur- 
face examined here. 



B. Nanocontact with adsorbed CO 

We will now describe the changes in the ballistic trans- 
mission induced by CO adsorption at the bridge site or 
at the atop site of the Au chain in the nanocontact. The 
transmissions are shown as a function of the scattering 
energy in Figs. [3ja-b) for the bridge geometry and in 
Figs.Efc-d) for the atop one (squares), together with the 
previously discussed transmissions of the clean nanocon- 
tact (diamonds) for the same values of d ss . 

When CO is at the bridge site, the conductance of 
the chain at low strains [Fig. ^l&)] is reduced to about 
0.63 Go, a value close to the fractional conductance peak 
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FIG. 4. (Color online) Transmission of a short chain with CO adsorbed at the bridge (left) or at the atop site (right): dependence 
on the smoothness of the surface/chain interface. The transmission of the (2v2 x 2\/2)R45 cell (abrupt interfaces, see Fig. [3} 
and that of the (3 x 3) cell (smoother interfaces, see insets) are compared for a selected value of d ss and are reported as a 
function of the scattering energy (squares and triangles, respectively). The central region of the smooth junction is built from 
the atomic positions optimized for the smaller cell (triangles) or by relaxing again the atoms in the constriction (circles). 



seen in Au nanocontacts experiments in presence of CO 
gas.— The transmission curve shows only a slight depen- 
dence on the scattering energy and decreases monoton- 
ically in the range of energies considered here. The 5d 
transmission peak seen below Ep in the pristine chain is 
suppressed here, but actually for energies slightly above 
or below it the transmission increases after CO adsorp- 
tion. The dependence of the transmission on the scat- 
tering energy around Ep is similar to that seen in Ref. 2 
for the tipless infinite chain geometries at medium/high 
strains, but the transmission values are lower by about 
10 % in the short chain geometries. At larger strains 
[Fig. [3jb)], the transmission of the bridge configuration 
is still characterized by a smoothly varying and mono- 
tonically decreasing behaviour as a function of energy, 
but shows slightly higher values with respect to the low 
strain configuration. The conductance grows to about 
0.72 Go and is therefore compatible with the conductance 
increase due to the contact stretching observed in the ex- 
perimental Au conductance traces at fractional conduc- 
tance values^ Previous conductance calculations for a 
3-atom chain between Au(lll) leads and a CO molecule 
at the bridge have found a larger conductance,— about 
0.9 Go, for a strain level which is roughly intermediate 
between the two considered here. This discrepancy may 
be ascribed to the shorter chain length^! or to the rather 
different functional used (B3LYP). 

When CO is adsorbed at the atop site, the conduc- 
tance reduction with respect to the clean nanocontact 
is much stronger than for the bridge site. In the low- 
strain configuration [Fig. [3jc)] , the transmission curve 
has a wide depression centered just above Ep and the 
conductance is about 0.08 Go, more than one order of 
magnitude smaller than in the clean nanocontact. Below 
Ep, the d peak is suppressed as in the bridge geometry, 
while at energies higher than 0.5 eV the s channel is less 
reflected and has a transmission which approaches that 
of the pristine chain. At higher strains [Fig. E^d)], the 



transmission dip shifts towards Ep causing a further low- 
ering of the conductance down to about 0.03 Go- There- 
fore, the strain dependence of the conductance in the 
bridge and in the atop configurations are completely dif- 
ferent, as already inferred from the transmission of the 
infinite chain geometries. 2 Moreover, the mechanism sup- 
pressing the transmission around the Fermi energy in the 
atop geometry is the same for both nanocontact and in- 
finite chain geometries. It can be regarded as a result of 
Fano-like destructive interference due to resonance scat- 
tering on the 5er a antibonding state brought up by CO 
adsorption and appearing right at Ep in the nanocontact 
geometry. A more complete picture of the connection be- 
tween the electronic structure features and the ballistic 
transmission of the CO/Au chain system was presented 
in a previous work using a tipless infinite chain model. 2 

These conclusions are not modified when considering 
"smoother" junctions with pyramidal tips on both sides 
of the chain. We use the previously relaxed "abrupt" ge- 
ometries to build a smoother junction by inserting the 
CO and the Au chain between the two terminal atomic 
planes of the tips (see insets of Fig. @|. These planes 
are made by four atoms each, arranged to form a square, 
and are spaced by d ss along z. In Fig. 01 we compare the 
transmission curves of the so-obtained "smooth" inter- 
face (triangles) with those of the corresponding "abrupt" 
interface (squares) for a selected value of d ss , both for 
the bridge and for the atop geometry. A quite good 
agreement can generally be observed below Ep, while 
above Ep there are larger discrepancies, especially for the 
bridge geometry. Nevertheless, the almost perfect corre- 
spondence of the conductance values confirms that the 
striking difference between the bridge and atop conduc- 
tances is not influenced by the precise atomic geometry 
chosen to model the wire/surface interface. We further 
optimized part of the structure by letting the 4 basal 
atoms on each side move along the longitudinal direc- 
tion (z axis) and completely relaxing the atomic posi- 
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FIG. 5. (Color online) Ballistic conductance as a function of 
the tips displacement: comparison between simulations (left 
panel) and experiment (right panel). Both LDA (open sym- 
bols) and GGA (filled symbols) conductances are shown as a 
function of d ss for the bridge (squares) and atop (circles) ge- 
ometries. d BB is shifted by a constant value d ss ,o to keep into 
account the experimental offset in the displacement measure- 
ment. The experimental conductance trace has been adapted 
from Ref. [y. The same scales for the conductance and the 
displacement are used in the two panels. 



tions of the chain and of the CO. The structural changes 
are rather small and the total energy is lowered by only 
0.3 eV or less. The transmissions obtained from these 
relaxed smooth interfaces (dots in Fig. |4|) do not differ 
appreciably from the non-relaxed structure in the bridge 
geometry, while the 5rj a transmission dip moves slightly 
closer to E-p in the atop geometry. In both cases, these 
conductance values are very similar to those presented 
above and do not affect the conclusions drawn from the 
analysis of the abrupt interfaces. 



V. DISCUSSION AND COMPARISON WITH 
EXPERIMENTAL DATA 

A reproducible behaviour in the conductance traces of 
Au nanocontacts in presence of CO was reported by a re- 
cent experiment^ a sharp reduction of the conductance, 
from the 1 Go value of the Au monatomic chain down 
to about 0.5 Go, followed by a slow increase upon fur- 
ther stretching of the nanocontact and by the final drop 
into the tunneling regime after a small elongation. It is 
reasonable to assume that, at some point of the pulling 
cycle, a CO molecule sticks to the Au chain giving rise 
to the sudden reduction of the conductance. Our theo- 
retical conductances for the bridge and atop geometries, 
together with their strain dependence, allow us to sim- 
ulate the conductance trace after CO adsorption at the 
bridge or at the atop site, and hence to make a direct 
comparison with the experimental traces. 

In Fig. 02 the conductances of the bridge and atop ge- 
ometries are shown as a function of the displacements 
between the Au surfaces and are juxtaposed with an ex- 
perimental conductance traced which presents the fea- 



tures described here above. We notice that the bridge 
configuration gives conductance values close to experi- 
mental ones and the correct dependence on strain, while 
the atop one gives the opposite behaviour with strain and 
too low conductance values at high strains. In the bridge 
geometry, the GGA and LDA conductances show a very 
good agreement with each other. In the atop geometry, 
the GGA conductance decreases more slowly than the 
LDA one, but the dependence on the contact stretching 
is the same. These conductance results, together with the 
strong energetic preference for the bridge site (Sec. IIIIj). 
are compatible with the presence of the fractional peak 
at about 0.6 Go in the Au conductance histograms and 
with the absence of a low conductance tail, which would 
be possible only in presence of an energetically favored 
atop geometry. 

However, we notice that the theoretical conductance of 
the bridge geometry reproduces the slope of the experi- 
mental conductance trace only qualitatively (Fig. [5]). We 
are not aware of other experiments reporting about the 
conductance slope of the CO/Au nanocontact system and 
we are not able to assess the degree of the experimental 
reproducibility, but we find anyway interesting to explore 
further this point. The discrepancy in the numerical val- 
ues of the slope might be related to the elastic response 
of the bulk leads to the external pulling force, which is 
present in the real nanocontact but has been neglected in 
the atomic relaxations of our model geometries. Indeed, 
a more realistic modeling of the structural modifications 
during the contact stretching process would require the 
relaxation of the atomic position in the atomic planes 
forming the leads, but this can be done only considering 
much larger supercells. 

It is possible to keep into account approximately the 
mechanical response of the leads to the external stress 
by treating them as ideal springs with a finite spring- 
constant k s M- We assume that the position of the elec- 
trodes is controlled at two opposite ends far away from 
the junction, separated by a distance L 3> d ss . The force 
balance between the ideal springs and the elastic response 
of the junction region can be expressed by the following 
equation*^ 



\UL 



d ss ) = 



dE tot (d S s) 
dd ss 



(1) 



where Et ot (d ss ) is the total energy of the relaxed 
nanocontact geometry with inter-surface distance d ss . 
By solving for L in Eq. ([1}, we can convert the inter- 
surface distance d ss into an equivalent tips displacement 
L(d S s', k s ), which includes the elongation of the leads and 
has only the electrode stiffness k s as external parameter. 
Notice that in the unrealistic assumption of infinitely stiff 
leads (k s = oo), the external stress is totally released in 
the nanocontact region that has been atomically relaxed 
(the short monatomic chain) and we have L = d ss . Ex- 
perimental estimates of k s in Au electrodes have been 
reported on the basis of the non-exponential dependence 
of the tunneling current on the distance: values in the 
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FIG. 6. (Color online) Top panel: ballistic conductance of the 
4Au-atom nanocontact with CO at the bridge site as a func- 
tion of the equivalent tips displacement L for different values 
of the electrode stiffness k s (in eV/A ) . The arrows point the 
start and end points indicated by the arrows in Fig. [5] and the 
solid line between those points approximates the experimen- 
tal conductance traced with a linearly increasing conductance. 
Bottom panel: total energy of the optimized nanocontact ge- 
ometry as a function of the inter-surface distance d ss . 



range 0.3 eV/ A to 3.7 eV/ A have been found, depend- 
ing on the junction realization^ 

In Fig. IS] (bottom panel), we report the total energy 
Etot(dss) of the relaxed LDA bridge geometry for some 
dss values in a range of interest and we fit this energy 
with a 4 th degree polynomial. The analytical deriva- 
tive of the energy fit is used to obtain L(d ss ;k s ) for a 
few values of the parameter k s and then express the (in- 
dependent conductance as a function of L (see Fig. [51 
upper panel). The case of infinitely stiff leads (k s = oo) 
corresponds to the conductance already shown in Fig. [5] 
and is reported here again for comparison. When taking 
a finite k s , a fraction of the term depending on the en- 
ergy derivative is added to d ss in i, thus the calculated 
conductance points are shifted to higher displacements 
and the theoretical conductance slope decreases. For k s 
values close to the experimental upper limit (4.0 eV/A ) 
the slope does not change significantly, but for k s val- 
ues between 0.5 eV/A and 1.0 eV/ A (a realistic range, 
according to the experiment) the calculated conduc- 
tance slope is in much better agreement with that of the 
experimental conductance traced (see thick solid line in 
Fig. [5]). This demonstrates that a more accurate descrip- 
tion of the leads is needed to obtain a correct slope of 
the conductance as a function of the contact stretch, but 
already good estimates can be obtained through this sim- 
ple elastic-response model. 



VI. CONCLUSIONS 



We have studied through density functional calcula- 
tions the adsorption of CO on Au monatomic chains in a 
model nanocontact geometry and its effects on the bal- 
listic conductance of the nanocontact as a function of the 
contact stretching. By comparing the adsorption energies 
of CO at the bridge and atop sites, we find that the bridge 
site is energetically favored at all levels of Au strain, as 
also found in a simpler model without bulk leads^ The 
chemisorption energies of CO on the short chain in the 
nanocontact geometry are comparable to those found in 
the infinite straight chains^ and a fair agreement in the 
variation of the energetics with strain is obtained when 
the displacement of the Au atom below the molecule in 
the infinite chain geometry is taken into account. This 
was not clearly predictable, because the finite length of 
the chain and the non-uniform Au-Au spacing are not 
encompassed by the infinite chain model, while in the 
nanocontact geometries these aspects are included and 
play a role in the structural modifications of the short 
chain after CO adsorption. For instance, when CO is at 
the bridge site, at low strains the chain bends towards the 
molecule, while at larger strains the Au-Au bond below 
CO elongates much more rapidly than the others. With 
CO atop, the chain forms a zigzag geometry, more pro- 
nounced at lower strains, where the Au atom right below 
the CO moves towards the molecule. 



The electron transmission across the Au nanocontact 
in presence of CO displays some important features which 
have already been found in the tipless geometries, most 
notably the transmission dip close to the Fermi level in 
the atop geometry which causes a strong suppression of 
the conductance. This dip is not present in the bridge 
geometry at those scattering energies, therefore the con- 
ductance reduction with respect to the pristine nanocon- 
tact is much smaller. Also the dependence of the bal- 
listic conductance on strain does not change with the 
inclusion of the tips, modeled here using either abrupt 
interfaces between the chain and the Au(001) surfaces or 
smoother pyramidal junctions: the bridge conductance 
is found to increase slightly with strain, compatible with 
the experimental findings, while the atop conductance 
drops rapidly down to zero as the contact is stretched 
because of the 5cr a transmission dip moving closer to E-p. 
With the inclusion of tips, the bridge geometry shows 
conductance values close to those of the experimentally 
observed structure which forms in the Au nanocontact 
after CO exposition. The slope of the experimental con- 
ductance with respect to the contact stretching can also 
be reproduced with reasonable accuracy by our calcu- 
lations if the tips displacement is computed taking into 
account the elastic response of the bulk leads through a 
simple model. 
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